Draft version October 1, 2008 

Preprint typeset using 1^1^^ style cmulatcapj v. 10/09/06 



THE UDF05 FOLLOW-UP OF THE HUDF: II. CONSTRAINTS ON REIONIZATION FROM Z-DROPOUT 

GALAXIES 1 

P. A. Oesch^, C. M. CAROLLO^ M. Stiavelli^, M. Trent:'', L. E. Bergeron^, A. M. Koekemoer^, R. A. Lucas^, C. 
M. Pavlovsky^, S. V. W. Beckwith'*, T. Dahlen'*, H. C. Ferguson^, Jonathan P. Gardner'', S. J. Lilly^, B. 

MOBASHER'\ N. PANAGIA^"''''' 
Draft version October 1, 2008 

ABSTRACT 

We detect three (plus one less certam) zgso-dropout sources in two separate fields (HUDF and 
NICP34) of our UDF05 HST NICMOS images. These z - 7 Lyman-Break Galaxy (LBG) candidates 
allow us to constrain the Luminosity Function (LF) of the star forming galaxy population at those 
epochs. By assuming a change in only Af, and adopting a linear evolution in redshift, anchored to 
the measured values at z ^ 6, the best-fit evolution coefficient is found to be 0.43 ± 0.19 mag per 
unit redshift (0.36 ± 0.18, if including all four candidates), which provides a value of M,(z = 7.2) = 
— 19.7 ± 0.3. This implies a drop of the luminosity density in LBGs by a factor of 2 — 2.5 over the 
~ 170 Myr that separate z ~ 6 and z ~ 7, and a steady evolution for the LBG LF out to z ~ 7, 
at the same rate that is observed throughout the z ~ 3 to 6 period. This puts a strong constraint 
on the star-formation histories of z ^ 6 galaxies, whose ensemble star-formation rate density must 
be lower by a factor 2 at ~170 Myr before the epoch at which they are observed. In particular, a 
large fraction of stars in the z ~ 6 LBG population must form at redshifts well above z ~ 7. The 
rate of ionizing photons produced by the LBG population decreases consistently with the decrease in 
the cosmic star formation rate density. Extrapolating this steady evolution of the LF out to higher 
redshifts, we estimate that galaxies would be able to rcionize the universe by z ~ 6, provided that 
the faint-end slope of the z > 7 LF steepens to a ~ —1.9, and that faint galaxies, with luminosities 
below the current detection limits, contribute a substantial fraction of the required ionizing photons. 
This scenario gives however an integrated optical depth to electron scattering that is ~ 2cr below 
the WMAP-5 measurement. Therefore, altogether, our results indicate that, should galaxies be the 
primary contributors to rcionization, cither the currently detected evolution of the galaxy population 
slows down at z > 7, or the LF evolution must be compensated by a decrease in metallicity and a 
corresponding increase in ionization efficiency at these early epochs. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: high-redshift — galaxies: 
luminosity function 



1. INTRODUCTION 

One of the most fundamental open issues in obser- 
vational cosmology is what sources reionized the early 
Universe. Since the z > 6 quasar Luminosity Func- 
tion (LF) is found to b e too shallow to provide enough 
ionizing photons [e.g. iFan et all 120011 : iMeiksinI l2005t 
IShankar fc Math"unl2007t ). the most likely culprits are 
thought to be galaxies. A key question is therefore, how 
the star forming galaxy population evolves beyond the 
putative end of rcionization at z ~ 6, i.e., well into the 
epoch of rcionization. 

Direct searches for Lya emission provide a powerful 
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way to study hig h-redshift galaxies fe.g. live et al.ll200"6t 
IStark et al.ll2007l : IOta eFallBoOS) : however, it is difficult 
to estimate the (continuum and) ionizing flux of Lya 
emitters, and thus their contribution to rcionization. A 
more fruitful approach is provided by the Lyman break 
technique. 

Thousands of Lyman-break galaxies (LBGs) have been 
identified between z ~ 3 and 6 from optical images, 
providing an estimate for the rest-frame UV LF at 
those redshifts (e.g. Oesch et al 2007, h ereafter O07; 
Bouwens et al 2007, hereafter B07 ; see alsolSteidel et al.l 
19991: iSawicki h Thompso^ 120061 : lYoshida et all l2006t 
Iwata et al.l l2007l ). The evolution of the LBG LF 
throughout this redshift range is however still being de- 
bated: especially at the bright end, the strong clustering 
of bright galaxies and the tiny area coverage of space- 
based images prevent setting strong constraints on the 
LF. Another important factor is that the number of spec- 
troscopically confirmed LBGs is still small due to the 
faintncss of the sources. Within these caveats, O07 and 
B07 have published independent, self-consistent analy- 
ses, which agree on a remarkably steady evolution of 
AA/* ~ 0.35 mag per unit redshift over the z = 3 — 6 
epoch. 

At redshifts z > 6, the rest-frame Lya line shifts into 
the near-infrared (NIR) and the relatively low efficiency 
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and small field of view of current NIR detectors have 
limited the identification of LBGs. Thus, in contrast to 
the large number of sources detected out to z 6, the 
numerous searches for z > 7 LBGs have led to only a 
handful of candidates (e.g. Bouwens et al 2008, hereafter 
B08: Mannucci et al.ll2007l:lB ouwens fc Illingworthl[200fi : 



TABLE 1 

Depth and Area of the UDF05 and HUDF NICMOS 
Observations 



lYan fc Windhorst I l2004HStanwav et al. 20081 Searches 
for intrinsically faint z>7 LBGs around lensing clusters 
also resulted only i n small numbers of candidates (e.g. 
Richard et all 120061 : iBradlev et all 12008; Ric hard et all 
2008|) , due to the small areas subject to the lensing mag- 
nification. 

One of the most comprehensive studies is the recent 
analysis of BOS, who have measured the z ~ 7 LBG 
LF from a large compilation of both HST and ground- 
based data, covering a total area of ~ 271 arcmin^ over 
the t wo fields of the GOODS survey (|Giavalisco et al.l 
[200l . A gain within the caveats imposed by probing, 
with small number statistics, only the bright-end of the 
LF, these authors report that the dimming of the LBG 
LF which has been measured out to z ^ 6 continues with 
no substantial change out to z ~ 7. 

A summary of the previous work in the context of 
reionization is that, albeit within some restrictive as- 
sumptions, the detected z ~ 6 LBG population appears 
to be just barely able to mainta in the Universe ion- 
ized (e .g. iFerguson et al.l |200^ iLehnert fc Bremeil 
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also the detected z ~ 7 LBGs fall short of being able 
to pr ovide the photons require d to reionize the Uni- 
verse (jBolton fc Haehneltl [2007[) . It seems that fainter 
galaxies must have played a very imp ortant role (e.g. 
lYan fc Windhorst l[200l iRichard et al.l[200l . 

Here we present our independent derivation of con- 
straints on the z ^ 7 LBG LF. This is desirable in order 
to independently check the different effects of various un- 
certainties in the process of estimating LFs, including the 
data reduction, object detection as well as estimation of 
effective selection volumina. We use our UDF05 data 
as well as the optical (Beckwith et al. 2007) and NIR 
(Thompson et al. 2005) data of the Hubble Ultra Deep 
Field (HUDF). The UDF05 survey is a 204-orbit HST 
Large Program of ultradeep ACS and NICMOS obser- 
vations of multiple fields, each located ~10' away from 
HUDF (see O07). The UDF05 was designed to image 
with ACS-WFC the two NICMOS parallel fields (here- 
after NICP12 and NICP34) that were acquired while the 
HUDF was observed with the ACS and thus facilitat- 
ing a reliable search for z > 6 galaxies in these fields to 
unprecedented depth. 

The combined NIR NICMOS dataset that we study in 
this paper consists of three fields, covering a total area 
of 7.9 arcmin^, with 5cr point-source magnitude limits 
varying between 27.9 - 28.8 in F160W (TJieo) and 27.6 - 
28.6 in FllOW (Jno). The 50% completeness limits vary 
between 26.9 and 28.1 mag (see Tab. [T]and i |3J|) . Each 
of the fields consists of one deep pointing and a slightly 
shallower flanking area (by ~ 0.6 mag). These data, also 
used in BOS, provide the deepest NIR data available to 
date. 

In !J2] we describe the data and the steps adopted to 
identify the z ~ 7 LBG candidates and constrain the 
LBG LF. In SJ3] we present our results and their impli- 
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within an aperture of 0" 3 radius'' 50% completeness limit, as 
estimated in section ^3.11 

cation for the evolution of the LBG LF into the reion- 
ization epoch, and in S|4] we discuss the contribution of 
the detected z ~ 7 galaxies to reionization. We adopt 
the concordance cosmology deflned by f^M = 0.3, f^A = 
0.7, 77o = 70 kms"^Mpc-\ i.e. h = 0.7. Where neces - 
sary we assume VLbh^ = 0.02265 (iHinshaw et al.l[200l 
Magn itudes are given in the AB system (|Oke fc GunnI 
IT983h . 



2. OBSERVATIONS AND SOURCE SELECTION 

2.1. The UDF05 and HUDF data 

We reprocessed the HUDF data using improved data 
reduction algorithms relative to the publicly released 
data, including an iterative sky subtraction scheme. All 
NICMOS data were drizzled to a 0'.'09 pixel scale. The 
optical ACS data were rebinned to this same scale and 
convolved with a Gaussian filter to match the width of 
the NICMOS point spread function (PSF). We used the 
i?i6o-band dat a to detect galaxie s with the SExtractor 
program (Berlin fc Arnouta Il996l ). and measured their 
colors in elliptical apertures scaled to 1.5 x their Kron 
radii, matched to the detection image. Total magni- 
tudes were measured in 2.5 Kron (AUTO) apertures and 
were corrected for flux losses in the wings of the PSF 
(0.07-0.10 mag). 

Only detections with signal-to-noise ratio (S/N) larger 
than 4.5 within 0'.'3 radius apertures in the iJigo-band 
were considered. Flux variations in apertures of different 
sizes, laid down on empty (sky) areas, were estimated in 
order to correct for pixel-to-pixel noise correlation in the 
images; the rms maps were then rescaled so as to match 
the noise in 0'.'3 radius apertures. Undetected fluxes were 
replaced with their 1—a upper limits. The zeropoints of 
the Jiio and i/ieo filters were corrected by 0.16 and 0.04 
mag, respectively, in order to account for non-linearity 
of the NICMOS NIC3 detector (see de Jong 2006). Note 
that any additional effect of detector non-linearity would 
be minimal as our sample covers only a small range in 
magnitudes. 

2.2. The z^cjQ-dropout candidates 

Candidate galaxies at redshifts between 6.S and S were 
selected from the SExtractor catalog according to the 
following color criteria: 

(^850 — Jiw) > 1-3 

(^850 - Jno) > 1-3 + 0.4(Jiio - Hiao) 
{Jiin — Hi6o) < 1-2 
S/N{Hieo) > 4.5 A S/N{Jno) > 2 
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Fig. 1. — The z — J vs. J — H color-color diagram, on which 
we highlight the selection criterion for z > 6.8 galaxies (gray 
shaded area). Evolutionary tracks of star-forming galaxies, ob- 
tained with Bruzual & Chariot (2003) models, are shown with solid 
lines in blue-to-green colors, corresponding to dust obscurations of 
i?(-B — y) = 0, 0.15, 0.3. Low redshift galaxy track s, up to redshift 
z = 4 , are obtained from the galaxy templates of IColeman et all 
II 19801) . and are shown as dash-dotted tracks in orange-to-red col- 
ors. Dwarf stars are expected to lie in the hatched region (Knapp 
et al. 2004, Burrows et al. 2006, Pickles 1998). The filled black 
circles correspond to the four z ^ 7 candidate galaxies identified 
in our UDF05+HUDF dataset. 

S/N{Vqoq) < 2 A S/Nims) < 2 

This selection is indicated as a gray shaded region in 
the z — J vs. J — H color-color diagram of Fig. [TJ 
The figure also shows tracks of different galaxy types 
as a function of redshift. A star- forming galaxy with 
metallicity 0.2Zq and a stellar age of 10^ yr is shown 
with three blue-to-green solid lines, corresponding to 
E{B — V) = 0, 0.15, 0.3 mag, respectively (using the ex- 
tinction law of Calzetti et al. 2000). Such star forming 
galaxies enter the selection window at ^ ~ 6.8. Lower 
redshift galaxies are shown with dash-dotted lines in 
orange-to-red colors. The hatched region in the diagram 
corresponds to the expected location of cool dwarf stars 
of type M to L (Knapp et al. 2004'^, Burrows et al. 2006, 
Pickles 1998). 

Note that the color selection that we have adopted 
in this study is the so- c alled "conservative" criterion in 
iBouwens fc lUingworthl (|2006f l. Relative to the less re- 
strictive color window that is used in BG8, our selection 
criterion slightly reduces the effective selection volume 
by excluding galaxies with redshifts between z ~ 6.5 and 
z ^ 6.8; however, it is more robust against low redshift 
interlopers, and in particular it excludes more efficiently 
a possible contamination by passively-evolving z 1.7 
galaxies with a pronounced 4000 A break. 

All SExtractor sources that satisfy the above color cri- 
terion were visually inspected; four z > 6.8 LBG can- 

^ We have used the theoretical models of Burrows et al. 2006 
to convert the photometry of Knapp et al. to the ACS/NICMOS 
filter set. 
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Fig. 2. — Images of our four 2s50-dropout candidates. The boxes 
are 3''5 in size. Left to right, the columns show images in Higg, 
JuOt 2g5o, «775, and Veoe- The sources are very compact. All but 
the faintest candidate (NICP34-80) are securely detected both in 
Hiao and Jiio. 



didates were identified after removal of spurious detec- 
tions (e.g., stellar diffraction spikes, edge artifacts, etc.). 
These four high-z galaxy candidates are shown as filled 
black circles in the color-color diagram of Fig. [TJ Images 
of these sources are shown in Fig. [21 their photometry 
is listed in Table [2l All four z-dropouts are rather faint 
(iJi6o ^ 27) and very compact (they are essentially un- 
resolved at the resolution of the NICMOS images). Two 
of the candidates lie in the shallow flanking region of 
the HUDF; the remaining two are located in the shal- 
lower area of the NICP34 data. As evident from Fig. 
[TJ three out of our four high-z candidates support rela- 
tively small amount of dust, E{B — V) — 0.15, while the 
last one (NICP 34-66) is consistent with slightly larger 
extinction. 

Three of our four sources are securely detected in both 
Jiio and HiQo- The faintest source (NICP34-80) has a 
S/N of 4.8 in the 0'.'3 iJieo aperture, and of only ~2 
in the Jno aperture. Thus NICP34-80 is a less secure 
high-z candidate. With the exclusion of NICP34-80, the 
remaining three objects are also found in the indepen- 
dent analysis of B08. One of the zgso-dropouts in our 
analysis (HUDF-48 0) had already been repo rted in previ- 
ous works feouwens et al.ll2004bl : IBouwens fc HHngworthl 
l2006t ICoe et ani2006D. and it has also been detected in 
Spitzer images ([Labbe et al.ll2006f) . with infrared colors 
in very good agreement with a star-forming galaxy at 
z ^ 7. The remaining Z85o-dropouts in our sample are 
either undetected (< 2a; NICP34-66/80)) or severely 
blended (HUDF-708) in the existing, publicly available 
Spitzer data (from the IRAC foUowup of GOODS for 
the HUDPS and from the SIMPLE survey for NICP349). 
The magnitude limits for the sources in NICP34 are quite 
bright and thus do not add any constraint on their pho- 
tometric redshifts. 

* http:/ /data. spitzer. caltech.edu/popular/goods/ 
^ http:/ /data. spitzer. caltech.edu/popular /simple/ 
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TABLE 2 

Photometry of the four, zsso-dropout sources 



ID 


a 


S 




S/N(Hi6o) 


■^110-^^160 


^SSO-'.'^llO 


FWHM 


HUDF-480 


3:32:38.81 


-27:47:07.2 


26.9 ±0.1 


7.2 


0.3 ±0.3 


2.4 ±0.7 


0'.'30 


HUDF-708 


3:32:44.02 


-27:47:27.3 


27.3 ± 0.2 


6.8 


0.5 ±0.3 


> 2.4 


0'.'29 


NICP34-66 


3:33:08.29 


-27:52:29.2 


26.9 ± 0.1 


12.5 


0.7 ±0.3 


> 1.6 


0'.'54 


NICP34-80'' 


3:33:10.63 


-27:52:31.0 


27.8 ±0.2 


4.8 


0.4 ±0.4 


> 1.8 


0'.'35 



This object could bo a spurious detection. We have detected it with a 2cr significance in Jiio, and it has not been confirmed by an 
independent reduction of the same data by R. Bouwens (2008, private communication). 



2.3. Sources of sample contamination 

There are five possible sources of contamination to the 
derived sample of z ~ 7 galaxy candidates; these are due 
to (i) spurious detections, (ii) cool dwarf stars, (Hi) low 
redshift galaxies with colors similar to z ~ 7 galaxies, {iv) 
high redshift supernovae, and (v) photometric scatter of 
low redshift galaxies into our selection window. In detail: 

{i) Except for the faintest of our sources, NICP34-80, 
all candidates are detected at more than 6 tr in Hiqq and 
> 3(7 in Jiio- It is therefore very unlikely that these 
high S/N sources are produced by peaks in the noise. 
In order to estimate the reality of NICP34-80, which is 
only detected at 2 cr in Jno, we repeated our analysis 
on 'negative images', obtained by multiplying the real 
images by —1 (after masking out all detected sources). 
One 'negative source' is detected with S/N(i?i6o)> 4.5, 
S/N(Jiio)> 2 and a pixel flux distribution similar to 
NICP34-80. Thus, we conclude that NICP34-80 is likely 
a spurious detection. 

(ii) Except for NICP34-66, whose profile has a full- 
width-half-maximum (FWHM) of 0'.'54, all sources are 
unresolved (FWHM - 0'.'29 - 0'.'35; see Table EJ. For 
comparison, s tars are measured to ha ve FWHM of 0'.'30 — 
0'.'39 (see also lThompson et al.ll2005l) . Thus, on the basis 
of their profiles, we cannot exclude the possibility that 
our candidates are stars. However, they all show colors 
redder than expected for ultracool dwarfs. Such stars 
should have at least weak detections in our ultra-deep 
Z850 data. Therefore it is unlikely that our candidates 
ar e dw arf stars. However, simulations of dwarf star LFs 
bv lBurgasse r (2004) predict that there will be ~ 0.5 stars 
of type L0-T8 in our survey in the magnitude range of 
Hab = 27 — 28. Thus, this is a non-negligible source of 
contamination. 

(Hi) Galaxies at z ~ 1.2 — 1.8 which emit a combina- 
tion of a red continuum (due to dust or old stellar age) 
and strong [OIII]A5007,A4959 and H/3 emission lines can 
exhibit colors which place them within our selection win- 
dow. The equivalent widths of the emission lines that 
are required are however rather extreme (> 500 A rest- 
frame), which makes this a rather implausible source of 
contamination. 

(iv) The 2.2 arcmin^ NICMOS observations for the 
UDF05 data were taken at a much earlier epoch than 
the corresponding ACS images; high redshift SNe could 
therefore also contaminate our sample. However, with a 
similar calculation as in BOS, we expect to find only 0.03 
SNe within the area of our NICMOS observations; this is 
thus a negligible source of contamination in our sample. 

(v) Photometric scatter could in principle cause some 
contamination as well. In particular, z ^ 1.7 galaxies 
with prominent 4000 A breaks, the closest to our selec- 



tion window, could be a significant source of contami- 
nation. In order to assess the importance of this effect, 
we run simulations in which we applied a suita ble photo- 
metric scatter to the HUDF source catalog of iCoe et alj 
(|2006f ). This catalog contains sources separately identi- 
fied both in the optical and in the infrared, and includes 
one of our z ^ 7 candidates. The catalog also contains 
photometric redshifts for all sources. We can thus study 
the redshift distribution of galaxies which enter our se- 
lection window after applying the photometric scatter. 
Certainly this test is limited by the reliability of the pho- 
tometric redshifts. However, it provides a worthy cross 
check on our candidates, and we found no low redshift 
source in any of the simulation runs. In general, requir- 
ing that an object is not detected in any of the optical 
bands efficiently reduces the contamination from lower 
redshift interloper galaxies. 

In summary, three of the four detected zsso-dropouts 
are very likely to be genuine z 7 LBGs. However, we 
stress that future observations will be needed to confirm 
these candidates as high-z sources. In the following we 
discuss the implications for galaxy formation and reion- 
ization of assuming that these three objects are truly 
galaxies at z 7; we also separately comment on the 
impact of including the fourth, less likely high-z candi- 
date in our z 7 LBG sample. 

3. RESULTS 
3.1. Expected vs. detected number of sources 

In order to estimate the completeness, C, and redshift 
selection probability, S, for our sample, we used the pro- 
cedure that we outlined in lOesch et al.l ()2007l ). i.e., we 
simulated data by inserting artificial galaxies into the 
real images, and run SExtractor on the simulated data 
with the same parameters adopted to extract the original 
catalogs. We adopted colors according to a Gaussian dis- 
tribution of continuum slopes that is measured for z ^ 6 
LBGs (/3 = -2.2 ± 0.2; Stanway et al. 2005), and a log- 
normal size distribution with a size scaling of (1 -I- z)~^ 
(iFerguson et al.ll2004 iBouwens et al.ll2004af) . 

Given a LF, the number of objects expected in a bin 
of a given magnitude, 0, is given by 

r r°° dV 

N{m)= dm' dz—(l){M[m',z])S(z,m')C(m') 

Jhin Jq dz 

where M[m,z] ^ m — K{z) — DM{z) is the absolute 
magnitude, m the observed magnitude, K{z) the K- 
correction term and DM the distance modulus. We 
adopted a 10^ yr old, 0.2Zq, star-forming template for 
the A'-correction from the observed i?i6o-band to rest- 
frame 2000 A, with E{B - V) = 15, consistent wit h 
best-fit SEDs for z ~ 5 LBGs (e.g. IVerma et al.llMl . 
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Note that correcting to 1400 A results in a dilFerence of 
only 0.07 mag. The total number of expected z ~ 7 
galaxies in our survey is then just the sum over all the 
fields, which all have different selection probabilities and 
completenesses due to the different depths and data qual- 
ity. 

Assuming no evolution in the LF from z 6 to 7, we 
would expect to detect 12 galaxies in our dataset (using 
the z ~ 6 parameters of B07); in contrast, only three, or 
at most four (if the least reliable candidate, N1CP34-80, 
is also included), are detected. In particular, we would 
expect about four galaxies brighter than Hiqq < 26.75. 
No such bright candidate is observed, which leads us to 
conclude, albeit with a limited significance, that there is 
evolution in the LBG LF from z ~ 6 to z ~ 7. From 
Poissonian statistics, the detection of four sources out of 
the expected 12 has a significance of 2.3 a; when cosmic 
variance is considered^°, the significance of the result is 
reduced down to 1.4 a. The same, i.e. a chance of 8%, is 
found from accurate beam tracing through a dark mat- 
ter simulation in which halos are populated with galaxies 
taking into account the specific geometry of our survey 
(jTrenti fc Stiavellil [20M l. Note however that, as com- 
mented above, the NICP34-80 source is likely a spurious 
detection; should this or any other high-z candidate in 
our sample turn out to be an artifact (or a lower redshift 
interloper), this would imply a stronger evolution of the 
LF over the z ~ 6 out to z ~ 7 period. The estimates 
above therefore provide a lower limit to the real evolution 
of the LF across the reionization boundary. 

3.2. The z ^ 7 LBG luminosity function 

The approach that we use to constrain the z ~ 7 LBG 
LF from our data is to fix the faint-end slope a and the 
normalization 0*, and search for evolution in M, only; 
this approach is motivated by the finding that M, is the 
only LF parameter which substantially varies over the 
z ~ 4 — 6 redshift range (B07). We adopt a linear evo- 
lution M^z) = M*(z = 5.9) +/3^{1 + z) and fit the only 
free parameter, to our four (three) z ^ 7 candidates. 
We use the parameters of B07 to anchor the evolution 
at redshift 5.9, i.e. we adopt (/),(z = 5.9) = 1.4 x 10"^ 
Mpc-^mag-i, a{z = 5.9) -1.74, and M,(z = 5.9) = 
-20.24 ±0.19. 

To determine /3,, we maximize the likelihood 
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where ^[a;. A] is the Poisson distribution with mean A 
evaluated at x\ Ndhsijni) is the observed number of ob- 
jects in the magnitude bin m^, and A'^'p are the expected 
number of objects for a given evolution parameter /3*. 

The best fit for /3* is 0.36 ± 0.18, which translates 
into M,(z = 7.2) = -19.77 ± 0.30. If we exclude the 
faintest, least reliable zgso-dropout candidate from the 
fit, the result changes only slightly, i.e., /3* — 0.43 ±0.19, 
or M,(z = 7.2) = -19.68 ± 0.31. Within the errors 
these values are identical, and we adopt a fiducial value 
of M,(z = 7.2) = -19.7 ± 0.3 at 2000 A. With this 
evolving LF, the redshift distribution of our sources is 

Our cosmic variance calculator is publicly available at 
|http://www. stsci.edu/~trenti/Cosmic Variance. html| 



10 




• this work 
n Bouwens et al. 2008 

* Mannucci et al. 2007 



-22 



-21.5 



-21 



-20.5 



-20 
Muv 



-19.5 



-19 



-18.5 



-18 



Fig. 3. — The z ~ 7.2 LBG LF derived from our four candidates 
(filled blue circles). The best-fit Schechter function, derived from 
maximizing the likelihood for the evolution parameter /3» , is shown 
as a dashed line. The dotted line is the same, but with a steeper 
faint end slope a = —1.9. Upper limits for empty magnitude are 
indicated with blue arrows. Open squares are from Bouwens et al. 
(2008), with whom we have three high-z candidates in co mmon. 
The upper Umit at M = -21.4 is from IMannucci et aT] (2007). 
The other lines show the LF at 2 ~ 5 (gray solid line; Oesch et al. 
2007), and 2 ~ 6 (black solid line; Bouwens et al. 2007). 



predicted to be z = 6.8 — 7.7 (within the 16th and 84th 
percentiles), with a mean z = 7.2. The resulting LF is 
shown in Fig. [3] as a dashed line; the figure assumes that 
all candidates lie at z ~ 7.2, with the effective volume 
estimated as V^?i{m) = dz^S{z,m)C{'m). 

Since the measured evolution is mostly driven by the 
non-detection of bright sources, the effect of varying the 
faint-end slope a is very small. Varying a by 0.16 (the 
formal l-cr error of the z ^ 6 LF) changes the value of 
(z = 7.2) by only ±0.03 mag. We have also tested how 
the LF would evolve if does not evolve while varies 
with redshift; this would result in a 0* value a factor 
-2.5 lower than at z - 6, i.e., 0* = (0.56 ± 0.50) x 10"^ 

B08 used a large sample of NIR data, covering an area 
of 261 arcmin^ (including the ultra-deep NICMOS data 
of this work) and found = —19.8 and (j)* = 1-1 x 10"'^ 
mag~"'^Mpc^^, in very good agr eement with our result 
(see also Fig. |3]). Furthermore, Mannuc ci et al.l ()2007l ) 
searched the ground-based VLT/ISAAC J, H, and Ks 
data on the GOODS field for zgso dropouts and, after 
removing probable dwarf stars, they found no candi- 
dates down to a UV rest-frame absolute magnitude of 
Afi5oo = —21.4. These authors inferred an upper limit 
on the LF down to this magnitude which we also report 
in Fig. [21 this upper limit is in good agreement with the 
extrapolation to brighter magnitudes of our z — 7 LF. 

4. IMPLICATIONS FOR GALAXY EVOLUTION 
AND REIONIZATION 

4.1. Evolution of the star formation rate density and 
star formation histories of z 6 LBGs 
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From the above LF we can compute the amount of evo- 
lution in the cosmic star formation rate density from z ^ 
6 out to z ~ 7. Integrating down to 0.2L,(z = 6), cor- 
responding to M = —18.5, we find a decrease in the lu- 
minosity density by 50%, i.e., log(p[erg/s/Hz/Mpc'^]) = 
25.5±0.2. Using the relati onship between UV continuum 
and star formation rate of iKennicuttI (|1998f ). this corre- 
sponds to a star formation rate density of io^2.32±o.23 
M© yr~^Mpc~'^. Assuming that the LF evolves only in 
density implies instead a reduction of 60%. Thus, inde- 
pendent from the specific assumption on the LF evolution 
(either in A/* or in (\)^\ down to this brightness limit the 
luminosity density, and thus the star formation rate den- 
sity, drops by a factor of ~ 2—2.5 over th e ~ 170 Myr that 
separa te the z = 6 and z = 1 epochs. iMannucci et al.l 
(|2007t l and BOS independently find similar results. Note 
that the conversion of Kennicutt is applicable for galaxies 
with continuous star-formation over timescales of a few 
100 Myrs and Salpeter initial mass function (IMF). They 
might differ for galaxies which arc experien cing stochas- 
tic bu rsts, as is assumed for high-z LBGs (jVerma et al.l 
l2007f ). However, larger differences are expected due to 
uncertainties in the IMF. 

The observed decline of the bright end of the LF from 
z ~ 6 out to z ~ 7 provides an upper boundary on 
the ensemble average star formation rate of z = 6 LBGs 
~ 170 Myr prior the epoch at which they are observed. 
Eyles et al. (2007) have investigated the ensemble aver- 
age star-formation history of i775-dropouts brighter than 
Z850 = 27 in the GOODS field, and suggested that this 
rises from z ~ 6 to a peak at z ~ 7.5, and declines beyond 
this redshift. Using Bruzual & Chariot (2 003) models 
and a ssuming the star formation history of lEvles et al.l 
(|200l . we find that the rest-frame 1500 A luminosity 
density should show an increase by a factor of ^ 2.5 from 
z ~ 6 up to 7, i.e., in striking contrast with the drop of 
the luminosity density that we report above. Dust cor- 
rections can not reconcile this disagreement, unless an 
evolution towards larger dust obscuration at higher red- 
shifts is invoked. 

This apparent contradiction could however be an arti- 
fact of the specific star formation history that has been 
adopted by Eyles at al. in the computation of the en- 
semble average. Generally, Spitzer data of i775-dropouts 
in the GOODS fields have revealed prominent 4000 A 
breaks in a large fraction of these sources (~ 40%), indi- 
cating u p to 90% of old stellar populations i n those z 6 
galaxies ()Yan et al.ll2005l : lEvles et al.ll2007| ). Our results 
on the z 7 LF provide global support to a scenario 
where a large fraction of stars in the z ^ 6 population 
form at redshifts well above z ~ 7. 

4.2. Do galaxies reionize our universe? 

Consistent with the decrease in cosmic star formation 
rate density, the rate of ionizing photons produced by 
L BGs must also decrease be tween z ^ 6 and ~ 7. 

iBolton fc Haehneltl (120071) pointed out tha t, with the 
estimated z ~ 7 LF of iBouwens et al.l ()2004bf ). the num- 
ber of ionizing photons falls short by almost an order 
of magnitude even to just keep the universe ionized at 
z ~ 7. Only models in which the star formation rate 
rises or stays constant beyond z ~ 6 are able to pro- 
vide enough ionizing photons. This result was however 
derived by including only galaxies with absolute magni- 
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Fig. 4. — Filling factor of ionized hydrogen as a function of 
redshift. The solid lines show the evolution of the filling factor cal- 
culated by assuming a faint end slope a = —1.7 A and an evolution 
of the galaxy LF by 0.35 mag per unit redshift, and by integrating 
down to different limiting SFRs (~ 1, 0.01, 0.0001 MQ/yr; lower 
gray, black, and upper gray lines, respectively). The dashed line 
assumes the same evolution in Af, but a = — 1.9. The dash-dotted 
line corresponds to no evolution of the LF from z ~ 6 out to z = 20; 
the dotted one assumes the same but with the z ~ 7 LF. If not 
stated otherwise all curves are integrated down to SFRs of 10"'^ 
MQ/yr. The corresponding electron scattering optical depths are 
also indicated. 



tudes brighter than M — —18 mag; this corresponds to 
a lower limit for the star formation rate of ~ 1 MQ/yr 
(see lower gray line in Fig. |4]). For LF faint-end slopes 
as steep as those that are measured, however, ^ 50% of 
the luminosity density comes from galaxies beyond the 
current detection limits at z '--^ 6 (Af > — 17 mag). Thus, 
photons from these very faint galaxies will add a sub- 
stantial contribution to the reionizing flux. 

Here we adopt the results of our independent analysis 
of the z ^ 7 LF, which, using substantially improved data 
products, confirms the drop in the star formation rate 
density from z = 6 to z = 7 suggested by Bouwens et al 
(2004b) . We follow the approach of Bolton & Hachnel: 
(|2007f ) to estimate the contribution of galaxies to reion 
ization. We include, however, the contribution of faint 



galaxies to reionization by decreasing the lower limit star 
formation rate to 10~^ MQ/yr in our calculations. 

Computations of the filling factor of ionized hydro- 
gen, QhiI) critically depend on the adopted values of 
the clumping factor C and escape fraction /esc. Sce- 
narios with larger clumping factors lead to more efficient 
recombinations and thus require more ionizing photons; 
lower escape fractions also require more photons. Both 
the clumping factor and the escape fraction arc however 
poorly known. To give an upper limit on the potential for 
reionization of the galaxy population, we adopt (possibly 
somewhat optimistic) val ues of C = 2, as found in simu- 
lations (see discussion in IBolton fc Haehn clt 2007), and 
/esc = 20%, the upper limit that has bee n inferred from 
both observations and simulations /e.g. [ Shaplev et al.l 
l20MISiaM"eraIll2007^ et al.ll2008bir We neglect 
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any possible mass dependence of /esc- 

Assuming (i) that the constant dimming of Af, with 
redshift continues beyond z ^ 7 and {ii) a faint-end slope 
a which remains constant at the value a = —1.74 that 
is measured at z ~ 6 (i.e., our "fiducial" evolution pa- 
rameters), the galaxy population would be able to ion- 
ize the universe only up to a fraction Qhii ~ 65% (see 
solid lines in Fig. |4]). The corresponding optical depth 
to electron scattering is r = 0.044, which is ~ 2.5a 
lower than the WM AP-5 value of r = 0.087 ± 0.017 
(jPunklev et al.| [2008[). Thus, this reionization scenario 
starts too late, ionizing too few atoms at high redshifts. 
Indeed the WMAP-5 data suggest that the Universe un- 
derwent an extended period of partial reionization, be- 
ginning as early as z ~ 20, simil ar to t he double reion- 
ization scenario proposed by Cei] (|2003f) . If by z ~ 9 the 
Universe is already substantially ionized by some other 
sources, the galaxy population would certainly be able 
to complete the reionization process by z ~ 6. 

Another possibility is that the z > 7 LP has a faint-end 
slope which is stee per than currentl y observed at 2 ^ 6 
(predicted e.g. bv lRvan et al. I [20071 ). The current esti- 
mate(s) do not allow one to rule out that the faint end of 
the LP at such early epochs is as steep as a = — 2, which 
would result in a diverging luminosity density. To assess 
whether the faint galaxy population can provide enough 
additional photons to complete reionization by z ~ 6, we 
adopt a value of a = —1.9 (see dotted line in Pig. [3]) and 
recompute the filling factor of neutral hydrogen under 
these modified assumptions. This scenario is shown in 
Pig. [Has a dashed line. Such a steep faint-end slope for 
the LP would result in an end of reionization by z 6, 
and in an integrated electron scattering optical depth of 
r = 0.049, which is still substantially lower than the cur- 
rent WMAP-5 value. Por comparison, we also show in 
Pig [H the evolution of the HII filling factor assuming no 



evolution in the LP from z ~ 6 (z ~ 7) to z = 20. In- 
tegrating the luminosity density again to 10^^ Mo/yr, 
this results in the completion of reionization by z 6.6 
(~ 5.4) and in an optical depth r = 0.083 (0.069). 

We therefore conclude that: (j) In order for the galaxy 
population to provide a substantial contribution to reion- 
ization, galaxies below current detection limits must play 
a significant role; (ii) Under this assumption, a galaxy 
population which evolves by 0.35 mag per redshift unit 
would be able to ionize the universe by z ~ 6, provided 
that the LP faint-end slope is steeper than the a = —1.74 
value that is measured at redshifts z < 6; {Hi) Even un- 
der this assumption, the resulting optical depth to elec- 
tron scattering is lower by ~ 2a relative to the WMAP- 
5 measurement, indicating that too few atoms are ion- 
ized at very high redshift. This finding provides evidence 
that, if reionization is primarily driven by galaxies, either 
the currently detected evolution of the galaxy population 
slows down at z > 7, or the LP evolution is compen- 
sated by a decrease in metallicit y and a corresponding in- 
crease in ionization efficiency dStiavelli et a l. 2004). Oth- 
erwise, other sources must substantially participate to 
reionize our Univer se like miniquasars or population III 
stars (e.g. jS huU & yenkatesanll2008t iMadau et ani2004l : 
ISokasian et al.ll200l l 
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